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ABSTRACT 

A complete sample of bright Swift Gamma-ray Bursts (GRBs) has been recently- 
selected by Salvaterra et al. (2011). The sample has a high level of completeness in 
rcdshift (90%). We derive here the intrinsic absorbing X-ray column densities of these 
GRBs making use of the Swift X-ray Telescope data. This distribution has a mean 
value of log(iV^f/cm -2 ) = 21.7 ± 0.5. This value is consistent with the distribution of 
the column densities derived from the total sample of GRBs with redshift. We find a 
mild increase of the intrinsic column density with redshift. This can be interpreted as 
due to the contribution of intervening systems along the line of sight. Making use of 
the spectral index connecting optical and X-ray fluxes at 11 hr {Pox), we investigate 
the relation of the intrinsic column density and the GRB 'darkness'. We find that 
there is a very tight correlation between dark GRBs and high X-ray column densities. 
This clearly indicates that the dark GRBs are formed in a metal-rich environment 
where dust must be present. 
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1 INTRODUCTION 

Long duration Gamma-ray bursts (GRBs) are thought to 
originate from the collapse of massive stars. Several lines of 
evidence points toward this conclusion, ranging from the as- 
sociation to type Ib/c supernovae (Woosley & Bloom 2006 
and references therein), to the occurrence of GRBs in the 
most luminous part of their host galaxies (Svensson et al. 
2010). The ambient medium in which GRBs explode is ex- 
pected to be denser than the interstellar medium and typical 
of star forming regions. The values of the absorbing column 
densities as measured in X-rays are high (Galama & Wijers 
2001; Stratta et al. 2004; Campana et al. 2006; Watson et 
al. 2007). An analysis of the intrinsic column densities of 
all Swift GRBs observed within 1,000 s and with a known 
redshift has been carried out by Campana et al. (2010). 
The selected sample consisting of 93 GRBs was biased. The 
distribution of the intrinsic X-ray absorption column den- 
sity is consistent with a lognormal distribution with mean 
log N H (z) = 21.9 ± 0.1 (90% confidence level). This distri- 
bution is in agreement with the expected distribution for 
GRBs occurring randomly in giant molecular clouds similar 
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to those within the Milky Way (Campana et al. 2006; Re- 
ichart & Price 2002). Looking at the distribution of X-ray 
column densities vs. redshift, there is a lack of non-absorbed 
GRBs at high redshift and a lack of heavily absorbed GRBs 
at low redshift. This might be the outcome of biases present 
in the sample. Looking at the distribution of X-ray col- 
umn densities versus redshift a lack of non-absorbed GRBs 
at high redshift and a lack of heavily absorbed GRBs at 
low redshift were found in previous studies (Campana et al. 
2010). The former might be explained in terms of more com- 
pact and dense star formation regions in the young Universe 
(or to a sizable contribution from intervening systems). The 
latter might be interpreted as due to a change in the dust 
properties with redshift, with GRBs at redshift z 2 — 3 
having a higher dust to gas ratio for the same X-ray col- 
umn density (e.g. different grain size or composition). This 
will naturally provide a lack of heavily (X-ray) absorbed 
GRBs at small redshifts. 

In the optical the presence of a large amount of ab- 
sorbing material is much less clear, since a large number of 
GRB afterglows are not affected by absorption (Kann, Klose 
& Zeh 2006; Schady et al. 2010; Zafar et al. 2011; but see 
Greiner et al. 2011 and Covino et al. 2011, in preparation). 
In this respect photoionisation of dust grains can play an 
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important role (Lazzati, Perna & Ghisellini 2001; Lazzati, 
Covino & Ghisellini; Draine & Hao 2002; Campana et al. 
2007). Moreover, the absorbing column densities measured 
in the optical based on damped Lyman-a absorption are 
a factor of ~ 10 lower than those measured in the X-ray 
band (Campana et al. 2010; Fynbo et al. 2009). This has 
been interpreted as due to photoionization of the surround- 
ing medium by GRB photons (Campana et al. 2006, 2007; 
Watson et al. 2007; Campana et al. 2010; Schady et al. 2011). 

The presence of a large amount of material is also testi- 
fied by the existence of 'dark' GRBs. There are several def- 
initions of dark GRBs. The easiest is that they do not show 
an optical counterpart (Fynbo et al. 2001). Since this defini- 
tion is clearly related to the sensitivity (and availability) of 
the instruments used for the follow-up a more general defini- 
tion is needed. Based on the predictions of the fireball model 
(Meszaros & Rees 1997) one can require that the optical to 
X-ray spectral index Box (i.e. the slope between the fluxes 
in the _R-band and at 3 keV at 11 hr after the burst) should 
be lower than 0.5 (Jakobsson et al. 2004). This will individu- 
ate optically sub-luminous bursts, i.e. fainter than expected 
from the fireball model. Alternatively, with the advent of 
Swift, X-ray spectral slopes were more easily available and 
a somewhat different definition was put forward by van der 
Horst et al. (2009) for which Box is shallower than fix —0.5. 

The darkness of a GRB can have different causes: it can 
be due to intrinsically optically faint GRBs, it can be due to 
absorption by intervening material within the host galaxy 
or it can be due to high redshift GRBs, thus being absorbed 
by the intergalactic medium. Several works have addressed 
this topic in the Swift era when a number of facilities al- 
lowed a quick follow-up of the afterglows. The fraction of 
dark bursts has been estimated to be ~ 25 — 50% according 
to Jakobsson's definition (Melandri et al. 2008; Roming et 
al. 2009; Cenko et al. 2009; Greiner et al. 2011; Melandri et 
al. 2011). It is now believed that the faint optical afterglow 
emission of dark bursts might be due to a moderate intrinsic 
extinction at moderate redshifts. Greiner et al. (2011) esti- 
mated a ~ 20% contribution from high redshift (z 4 — 5) 
GRBs to the dark population only. 

Salvaterra et al. (2011, see also Nava et al. 2011) se- 
lected a complete sample of bright GRBs based on optical 
observability (Jakobsson et al. 2006) and Swift BAT peak 
flux P 2.6 ph s" 1 cm -2 . The sample consists of 58 GRBs 
and it is complete in spectroscopic redshift at 90% (with 95% 
of GRBs having some constraints on the redshift) . This sam- 
ple offers the opportunity to study in an unbiased way the 
distribution of the X-ray column densities and its relation 
to GRB darkness. 

The paper is organised as follows. In section 2 we derive 
the X-ray absorbing column densities for the Salvaterra's 
sample and briefly describe how the slope Box has been 
computed for each GRB of the sample. In section 3 we dis- 
cuss our findings and in section 4 we draw our conclusions. 



2 X-RAY ABSORBING COLUMN DENSITIES 
AND SPECTRAL SLOPE BETWEEN 
OPTICAL AND X-RAY BANDS 

The intrinsic column densities were computed using the au- 
tomated data products provided by the Swift/XRT GRB 
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Figure 1. Column density (Nj{) as derived from the X-ray data 
of the GRBs in our complete sample. The dotted (orange) his- 
togram shows the distribution for GRBs with z < 1.7 and the 
dashed (light blue) histogram shows the distribution for GRBs 
with z > 1.7. 



spectra repository (Evans et al. 2009). The archive has been 
recently updated by reprocessing all the on-line GRB data 
products using the latest software and calibration (Evans 
2011). Therefore some of these values overwrite those re- 
ported in Campana et al. (2010). In Table 1 we list the col- 
umn density value at the host galaxy redshift Nh{z). These 
are obtained fitting an absorbed power law model to the 
data in the specified time interval when there are no strong 
spectral variations. The absorption component is modeled 
with PHABS within XSPEC. We consider two components one 
Galactic (held fixed) and a component at the redshift of the 
GRB (ZPHABS). The Galactic column density for each burst 
is provided by the Leiden/ Argentine/Bonn (LAB) Survey of 
Galactic HI (Kalberla et al. 2005). For those GRBs without 
redshift, we fix the redshift of the free absorption component 
to zero, so that the resulting value provides a lower limit to 
the intrinsic column density. 

In the next sections we will also use the Box index. 
This index is computed as the spectral index connecting the 
_R-band flux and the (unabsorbed) 3 keV flux at 11 hr from 
the trigger. The collection of indexes and limits for the burst 
in our sample can be found in Melandri et al. (2011). 



3 DISCUSSION 

3.1 Distribution of absorbing column densities 

The distribution of the X-ray column densities is shown in 
Fig. 1. The overall distribution can be well described by a 
Gaussian with (logarithmic) mean 21.7 and standard devia- 
tion 0.5 (the median of the distribution is very close to the 
mean). This is consistent with the column density distri- 
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Table 1. Column densities for bright Swift GRBs in the Salvaterra et al. (2011) sample. 



GRB z N H (z) T iV H (Gal) Time interv. Comments 

(10 21 cm- 2 ) (10 20 cm" 2 ) (s) (exp. time ks) 



050318 


1.44 


7 +0.6 
U - ' -0.5 


1.97 ±0.07 


1.9 


3000 - 7 x 10 4 


PC (23.5) 


050401 


2.90 


18.2+ 2 j 


1.84 ±0.05 


4.4 


200 - 2000 


WT (1.4) 


050416A 


0.65 


6 4+ ' 9 

° -0.5 


2.12 ±0.06 


2.4 


400 - 3 x 10 5 


PC (92.3) 


050525A 


0.61 


2-0±g-g 


2.12 ±0.16 


9.1 


6000 - 3 x 10 5 


PC (18.7) 


050802 


1.71 




1.86 ±0.12 


1.9 


400 - 2000 


PC (1.5) 


050922C 


2.20 


q c+2.2 
°-° — 2 1 


2.24 ±0.11 


5.4 


400 - 3 x 10 5 


PC (7.3) 


060206 


4.05 


14 2 +16 '° 
± ^* z — 9 8 


2.04 ±0.28 


0.9 


700 - 10 5 


PC (4.5) 


060210 


3.91 


24 -6 + 35 


2.09 ±0.05 


6.1 


3000 - 10 6 


PC (76.8) 


060306 


3.50 




1.86 ±0.26 


3.4 


300 - 925 


PC (0.6) 


060614 


0.13 


0.331°^ 


1.90 ±0.06 


1.9 


4400 - 10 6 


PC (28.4) 


060814 


1.92 


20.9+1J 


1.94 ±0.07 


2.3 


250 - 500 


WT (0.2) 


060904A 


- 


> 2.07 


3.53 ±0.43 


1.3 


185 - 225 


WT (0.1) 


060908 


1.88 


6 2+ 2 ' 8 


2.12 ±0.18 


2.3 


200 - 10 s 


PC (11.6) 


060912A 


0.94 


o n+1.5 
a - z -1.3 


2.02 ±0.18 


3.9 


200 - 2000 


PC (1.7) 


060927 


5.47 


< 36 


1.94 ±0.16 


4.6 


100 - 10 4 


PC (3.5) 


061007 


1.26 


r- 1+0.3 
°- i -0.3 


1.85 ±0.02 


1.8 


90 - 2000 


WT (1.9) 


061021 


0.35 


0.73t°' 2 


1.99 ±0.03 


4.5 


3000 - 3 x 10 5 


PC (83.1) 


061121 


1.31 


5 4+ ' 8 
°'^-0.5 


1.88 ±0.05 


4.0 


600 - 3 x 10 5 


PC (43.1) 


06 1222 A 


2.09 


44.811;* 


2.11 ±0.06 


9.0 


3 x 10 4 - 2 x 10 5 


PC (49.5) 


070306 


1.50 


26.8115 


1.88 ±0.12 


2.9 


10 4 - 4 x 10 4 


PC (6.0) 


070328 


- 


2 4+ ' 2 


2.24 ±0.05 


2.6 


350 - 1000 


WT (0.7) 


070521 


1.35 


5411? 


1.78 ±0.20 


2.9 


3000 - 10 4 


PC (1.9) 


071020 


2.15 


6.811-i 


1.87±0.07 


5.1 


70 - 300 


WT (0.2) 


071112C 


0.82 


I 4+0-5 


1.82 ±0.08 


7.4 


90 - 280 


WT (0.2) 


071117 


1.33 


10.9 +2 ' l 


2.05 ±0.18 


2.3 


2900 - 6.2 x 10 4 


PC (19.0) 


080319B 


0.94 


1 7+0. i 
1 - ' -0.1 


1.78 ±0.02 


1.1 


800 - 2000 


WT (1.7) 


080319C 


1.95 


5.5_ 2 3 


1.61 ±0.10 


2.2 


200 - 3 x 10 5 


PC (4.1) 


080413B 


1.10 


1 O+0.6 


1.97 ±0.07 


3.1 


6 x 10 3 - 10 6 


PC (40.6) 


080430 


0.77 


O K+0.8 

°- o -0.6 


2.03 ±0.10 


1.0 


5000 - 3 x 10 5 


PC (10.9) 


080602A 


1.40 


6 7+ 2 ' 4 

—2.1 


2.01 ±0.15 


3.5 


200 - 800 


PC (0.6) 


080603B 


2.69 


7 q+2.9 
'■ -2.7 


1.84 ±0.10 


1.2 


100 - 250 


WT (0.2) 


080605 


1.64 


9 0+ ' 9 


1.76 ±0.04 


6.7 


100 - 800 


WT (0.6) 


080607 


3.04 


22.8l4;2 


2.14 ±0.10 


1.7 


4000 - 6 x 10 4 


PC (19.7) 


080613B 


- 


> 0.5 


1.31 ±0.12 


3.2 


105 - 190 


WT (0.1) 


080721 


2.59 


10 4+ 06 


1.81 ±0.02 


6.9 


100 - 2000 


WT (1.3) 


080804 


2.20 


1 4+1-9 
— 1 1 


1.82 ±0.09 


1.6 


200 - 3 x 10 5 


PC (12.6) 


080916A 


0.69 


S.Olfj 


2.26 ±0.15 


1.8 


2 x 10 4 - 10 6 


PC (172.9) 


081007 


0.53 


4.81H 


2.04 ±0.18 


1.4 


5000 - 4 x 10 5 


PC (9.7) 


081121 


2.51 


1-9111 


1.93 ±0.06 


4.0 


3000 - 2 x 10 6 


PC (36.3) 


08 1203 A 


2.05 


45 +i.i 


1.74 ±0.05 


1.7 


200 - 600 


WT (0.4) 


081221 


2.26 


26.1 + |-| 


2.04 ±0.09 


2.0 


300 - 500 


WT (0.2) 


081222 


2.77 


6.011-i 


1.96 ±0.04 


2.2 


60 - 1000 


WT (0.8) 


090102 


1.55 




1.73 ±0.13 


4.1 


700 - 7 x 10 4 


PC (1.8) 


090201 


< 4 


> 3.85 


2.01 ±0.16 


4.9 


3000 - 6000 


PC (1.7) 


090424 


0.54 


4.i±S-.g 


1.94 ±0.08 


1.9 


2000 - 3 x 10 6 


PC (14.9) 


090709A 


< 3.5 


> 1.82 


2.04 ±0.13 


6.6 


4000 - 1.5 x 10 4 


PC (4.4) 


09U715B 


3.00 


7 - 6 -2.8 


2.01 ± 0.10 


1.3 


A i \: \i \ i r\5 

40U0 — 10 


FC (28.8) 


090812 


2.45 


1 2 0+ 7 ' 2 


2.10 ±0.23 


2.3 


10 4 - 7 X 10 4 


PC (7.7) 


090926B 


1.24 


13.9+H 


1.97 ±0.08 


1.9 


130 - 300 


WT (0.2) 


091018 


0.97 




1.84 ±0.12 


2.8 


150 - 1000 


PC (0.9) 


091020 


1.71 


5-811:5 


1.82 ±0.10 


1.4 


200 - 400 


WT (0.2) 


091127 


0.49 


0.76+g-f 


1.80 ±0.11 


2.8 


6000 - 2 x 10 4 


PC (2.0) 


091208B 


1.06 


s q+ 4 - 3 

8 - d -3.4 


2.16 ±0.27 


4.9 


200 - 600 


PC (0.4) 


100615A 




> 10.1 


2.43 ± 0.32 


3.3 


200 - 2000 


PC (1.4) 


100621A 


0.54 


18.0+1' 2 


2.86 ±0.09 


2.9 


130 - 200 


WT (0.1) 


100728B 


2.11 


4-3lli 


2.08 ±0.18 


6.2 


4000 - 3 x 10 4 


PC (8.5) 


110205A 


2.22 


3.5111 


2.11 ±0.09 


1.6 


5000 - 5 x 10 4 


PC (18.8) 


110503A 


1.61 


q r-q+0.67 


1.80 ±0.05 


2.6 


200 - 700 


WT (0.5) 



Errors and upper limits are at 90% confidence level obtained with a A\ 2 = 2.71. 
© 0000 RAS, MNRAS 000, 000-000 



4 S. Campana et al. 




Figure 2. Column density distribution versus redshift. Error bars 
have been computed from X-ray spectral fitting and are at 90% 
confidence level. Upper limits are at 90% confidence level. Light 
grey lines represent the possible column density values for GRBs 
without redshift, scaling the z = value with (1 + z) 2,6 within 
the range of allowed redshifts (Galama & Wijers 2001). Darker 
(orange) region limited by a continuous line marks the mean con- 
tribution in the observed column density jVjy(z) resulting from 
intervening systems as a function of redshift. The lighter (light 
orange) region limited by a dotted line marks the maximal 90% 
line of sight contribution of intervening systems to the observed 
column density. The dashed line marks the mean contribution to 
Njj (z) in the case of doubling the population of intervening sys- 
tems as suggested by the comparison of quasars and GRB studies 
of intervening systems (see text, section 3.1). 



bution obtained considering all GRBs with known redshift 
(Campana et al. 2010). The distribution of column densities 
as a function of redshift is shown in Fig. 2. It is apparent 
from Fig. 2 that there is a trend of increasing column den- 
sities with redshift. However, the lack of mildly absorbed 
GRBs at high redshift is not statistically overwhelming. To 
see if there is a real difference we cut the sample at the mean 
redshift (z = 1.7), and make a Kolmogorov-Smirnov (KS) 
test on the two distributions. We obtain a probability of 9% 
that the two distributions come from the same parent popu- 
lation. No firm conclusions can therefore be drawn. A cut at 
z = 1 results in a 0.5% KS probability. This might indicate 
a difference in the intrinsic absorption column densities (see 
below). With respect to a non-complete sample (Campana et 
al. 2010), we note that the high-column density region at low 
redshift is here more populated. This indicates a bias present 
in the non-complete sample: it is difficult to obtain the red- 
shift of highly-absorbed low-redshift GRBs (likely due to a 
higher optical absorption). 

The lack of mildly-absorbed high redshift GRBs re- 
mains. This has been interpreted as due to the presence 
of absorbing matter along the line of sight not related to 
the GRB host galaxy. This can be either diffuse (i.e. lo- 



cated in diffuse structures like the filaments of the Warm- 
Hot Intergalactic medium - WHIM, Behar et al. 2011) 
or concentrated into intervening systems (i.e. galaxies or 
clouds along the line of sight, Campana et al. 2010). Based 
on quasar studies (Wolfe et al. 2005; Peroux et al. 2003; 
Noterdaeme et al. 2009; Prochaska, O'Meare & Worseck 
2010), we can evaluate the contribution of the interven- 
ing systems to the observed column density by simulat- 
ing their distribution. In particular, we assumed a num- 
ber distribution of intervening systems, based on damped 
(and sub-damped) Lyman-a systems, with a redshift de- 
pendence n(z) oc (1 + z) ' 26 for z 2.3 and n(z) tx 
(1 + z) 2,78 for 2 > 2.3. For the contribution of the inter- 
vening systems in terms of absorption we adopted a three- 
fold equation (measuring column densities in cm -2 units): 
for logiVif < 18.2 we consider n(logiVff) oc logiVn -1 ' 9 , for 
18.2 < logTVj/ < 20.6 n(logiV H ) oc logA^ -08 , and for 
logiVff > 20.6 n(logTVff) oc logiNfe -1 ' 4 (Wolfe et al. 2005; 
Peroux et al. 2003). The contribution of each intervening 
system is weighted considering the system as if it was at 
the redshift of the GRB, i.e. weighting its intrinsic column 
density as ((1 + zgrb)/(1 + z DL a)) 2S (where zgrb is the 
redshift of the GRB and zdla the redshift of the interven- 
ing system) . We set up a MonteCarlo simulation considering 
systems in the log(jVjj /cm" 2 ) = 17.2 — 22 range and prob- 
ing 10,000 lines-of-sight. Because the GRB column densities 
were calculated assuming solar metallicity, for comparison 
we also assumed solar metallicities for the evaluation of the 
contribution of the intervening systems. It is important to 
note that GRB hosts have typically sub-solar metallicities, 
in which case the assumption of solar metallicity would lead 
to the equivalent column density being underestimated. All 
data points in Fig. 2 can therefore be (usually) considered 
as being lower-limits on Nh- The 90% confidence envelope 
of the simulated line-of-sights as a function of redshift is 
shown in Fig. 2 (dashed line). The average contribution is 
also shown (continuous line and dark orange region). These 
calculations clearly provide just an indicative estimate and 
are subject to uncertainties related to modeling the number 
density evolution in redshift of these systems (e.g. Ribaudo, 
Lehner & Howk 2011). It is apparent from Fig. 2 that a GRB 
lying along the 'average' line of sight experiences a too low 
increase of the intrinsic column density due to intervening 
systems with respect to the observed increase of GRB col- 
umn densities with redshift (even if a less favorable line of 
sight might fully account for the observed increase at high 
redshifts) . 

Studies on strong intervening systems in quasars and 
GRB spectra have shown a larger occurrance of intervening 
systems in the latter lines of sight (e.g. Prochter et al. 2006). 
These systems are mainly identified through strong, rest- 
frame equivalent width (EW) > 1 A Mg II AA 2796, 2803 
absorption lines. In a study with a large sample of GRBs, 
Vergani et al. (2009) confirm the presence of this effect and 
set the discrepancy to a factor of ~ 2. Even if the reason for 
this discrepancy is still not fully understood, Budezynski & 
Hewett (2011) show that this discrepancy is likely related to 
a lack of quasars heavily absorbed along their line of sight. 
Given this observational result, we artificially increased the 
number of intervening systems based on quasar studies by 
a factor of 2. The resulting mean contribution derived from 
the intervening systems is shown in Fig. 2 with a dashed line. 
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Figure 3. Column density (Nh) as a function of the spectral 
index Pox (Melandri et al. 2011) for the GRBs in our complete 
sample. Limits on the column density and Pox arc indicated with 
arrows. The dashed line for Pox = 0.5 divides dark and non-dark 
GRBs according to Jakobsson et al. (2004). Open squares (filled 
circles) indicate dark (non-dark) GRBs according to van der Horst 
et al. (2009). GRBs with lower limits on the column density (i.e. 
without rcdshift information) are not shown. 



This line follows nicely the increase of the intrinsic column 
density with redshift providing a plausible explanation for 
this effect. 



3.2 Relation between X— ray absorption and GRB 
'darkness' 

Given our complete sample of bright GRBs we investigate 
the connection between the X-ray absorption and the GRB 
darkness. The GRB darkness can be caused by several ef- 
fects that can be divided into two main classes: i) intrinsic, 
i.e. due to some physical mechanism hampering the optical 
emission or ii) environmental, i.e. due to intrinsic absorption 
and/or to high redshift. Considering the Pox values com- 
puted in Melandri et al. (2011), there are 19 GRBs in our 
complete sample that can be classified as dark according to 
Jakobsson et al. (2004) and 12 according to van der Horst 
et al. (2009). Out of them, 4 (common to both definitions) 
do not have any redshift information. In Fig. 3 we show the 
distribution of Pox as a function of the intrinsic column 
density for the GRBs of our complete sample (we did not 
include the few GRBs without a redshift determination). It 
is apparent that for bursts with Pox < 0.5 (i.e. dark ac- 
cording to the Jakobsson's definition) all but three (all with 
0.45 < Pox = 0.5) have an intrinsic column density larger 
than \og(NH /cm -2 ) £ 22. This same Nh limit is valid for 
all the bursts that are dark according to the van der Horst's 
definition (Fig. 3). 



Comparing the intrinsic column density distribution of 
dark and non-dark GRBs (taking only the ones with known 
redshift, i.e. 15-38 and 8-45 for the Jakobsson's and van der 
Horst's definition, respectively) we obtain a KS probability 
of 2 x 10~ 6 (4.8 a) for the Jakobsson's definition and 1 x 10~ 5 
(4.4 a) for the van der Horst's definition (the lower value is 
due to the smaller number of dark GRBs). We also note 
that if the 4 GRBs classified as dark and without redshift 
information would have a redshift in line with the mean 
of the sample, then they would have an intrinsic column 
density Nh(z) XL 10 22 cm~ 2 . These results indicate that the 
intrinsic absorption as evaluated in the X-ray band is highly 
correlated with the darkness of a GRB. 



4 CONCLUSIONS 

Salvaterra et al. (2011) selected a complete sample of bright 
GRBs with a high degree of completeness in redshift. In a 
series of papers we investigate the impact of this sample 
on GRB studies. Here we focus on the properties of the 
sample with respect to the intrinsic X-ray absorption. We 
found that the intrinsic column density distribution of our 
complete sample is consistent with the total distribution of 
column densities presented in Campana et al. (2010). The 
mean of the two distributions are in fact 21.7 ± 0.5 and 
21.9 ± 0.1, respectively. This likely indicates that the GRB 
brightness, as well as any other bias present in the total 
sample of GRBs with redshift (e.g. dust), does not heavily 
influence the total distribution of intrinsic column densities. 

At variance with the total distribution presented in 
Campana et al. (2010), we see in the complete sample pre- 
sented here that the region at high column densities and low 
redshift is now more populated by GRBs. This clearly re- 
veals a bias present in the non-complete sample, where this 
region is heavily underpopulated due to the lack of a redshift 
determination of dark bursts. 

Even if not statistically compelling there is an increase 
of the intrinsic column density with redshift (this is more ap- 
parent in the full sample of GRBs with redshift, Campana 
et al. 2010). We evaluate the mean contribution to Nh(z) 
due to the intervening systems along the GRB line of sight. 
We find that, if we take into account the larger number of 
observed systems affecting the line of sight of GRBs with 
respect to the quasar one (Vergani et al. 2009), the popu- 
lation of sub-Damped Lyman-a and Damped Lyman-a sys- 
tems can account for the increase with redshift of Nh(z). 
It would be interesting to confirm this directly through the 
study of high-0 GRB lines of sight. Unfortunately this effect 
plays a significant role at very high redshift, where the num- 
ber of GRB afterglow spectra is very low. It is indeed difficult 
to measure absorption from Lyman-a. However, the column 
density of neutral gas can still be traced by weakly ionised 
metal lines (e.g. Zn II, Si II), which in fact is a more logical 
method of comparing absorption in X-rays and the opti- 
cal, given that the X-rays are absorbed by metals and not 
neutral hydrogen (e.g. Schady et al. 2011). The X-shooter 
instrument mounted at the ESO/VLT offers the best oppor- 
tunities for these studies. 

Making use of the Pox computed by Melandri et al. 
(2011), we found a strong correlation between GRB dark- 
ness and X-ray absorbing column densities. Since metals are 



© 0000 RAS, MNRAS 000, 000-000 



6 S. Campana et al. 



Wolfe, A. M., Gawiser, E., Prochaska, J. A. 2005, ARA&A, 43, 
861 

Woosley, S., Bloom, J. S. 2006, ARA&A, 44, 507 
Zafar, T., Watson, D., Fynbo, J. P. U., Malesani, D., Jakobsson, 
P., de Ugarte Postigo, A. 2011, A&A, 532, A143 



5 ACKNOWLEDGMENTS 

SC thanks Darach Watson and Phil Evans for useful con- 
versations. This work has been supported by ASI grant 
1/004/11/0. This work made use of data supplied by the UK 
Swift Science Data Centre at the University of Leicester. 



REFERENCES 

Behar, E., Dado, S., Dar, A., Laor, A. 2011, ApJ, 734, 26 
Budzynski, J. M., Hewett, P. C. 2011, MNRAS, 416, 1871 
Campana, S., Thone, C. C, de Ugarte Postigo, A., Tagliaferri, 

G., Moretti, A., Covino, S. 2010, MNRAS, 402, 2429 
Campana, S., et al. 2006, A&A, 454, 113 
Campana, S., et al. 2007, ApJ, 654 L17 
Cenko, S. B., et al. 2009, ApJ, 693, 1484 
Draine B. T., 2003, ARA&A, 41, 241 
Draine, B. T., Hao, L. 2002, ApJ, 569, 780 
Evans, P. A., et al. 2009, MNRAS, 397, 1177 
Evans, P. A. 2011, GCN 12250 
Fynbo, J. P. U., et al. 2001, A&A, 369, 373 
Fynbo, J. P. U., et al. 2009, ApJS, 185, 526 
Galama, T. J., Wijers, R. A. M. J. 2001, ApJ, 549, L209 
Greiner, J., et al. 2011, A&A, 526, A30 
Jakobsson, P., et al. 2004, A&A, 427, 785 
Jakobsson, P., et al. 2006, A&A, 447, 897 

Kalberla, P. M. W., Burton, W. B., Hartmann, D., Arnal, E. M., 
Bajaja, E., Morras, R., Poppel, W. G. L. 2005, A&A, 440, 
775 

Kann, D. A., Klose, S., Zeh, A. 2006, ApJ, 641, 993 
Lazzati, D., Covino, S., Ghisellini, G. 2003, MNRAS, 330, 583 
Lazzati, D., Perna, R., Ghisellini, G. 2001, MNRAS, 325, L19 
Melandri, A., et al. 2008, ApJ, 686, 1209 

Melandri, A., et al. 2011, MNRAS in press (paper IV, 

larXiv:1112.4480H 
eszaros, P., Rees, M. J. 1997, ApJ, 476, 232 

Nava, L., et al. 2011, MNRAS in press (paper II. IarXiv:1112.4470H 
Noterdaeme, P., Petitjean, P., Ledoux, C, Srianand, R. 2009, 

A&A, 505, 1087 
Peroux, C, et al. 2003, MNRAS, 345, 480 
Prochter, G. E., et al. 2006, ApJ, 648, L93 

Prochaska, J. X., O'Meara, J. M., Worseck, G. 2010, ApJ, 718, 
392 

Reichart, D. E., Price, P. A. 2002, ApJ, 565, 174 

Ribaudo, J., Lehner, N., Howk, J. C. 2011, ApJ, 736, 42 

Roming, P. W. A., et al. 2009, ApJ, 690, 163 

Salvaterra, R., et al., 2011, ApJ in press (paper I, 

larXiv:1112.1700H 
Schady, P., et al. 2010, MNRAS, 401, 2773 

Schady, P., Savaglio, S., Kruhler, T., Greiner, J., Rau, A. 2011, 

A&A, 525, A113 
Stratta, G., Fiore, F., Antonelli, L. A., Piro, L., De Pasquale, M. 

2004, ApJ, 608, 846 
Svensson, K. M., Levan, A. J., Tanvir, N. R., Fruchter, A. S., 

Strolger, L.-G. 2010, MNRAS, 405, 57 
van der Horst, A. J., et al. 2009, ApJ, 699, 1087 
Vergani, S. D., Petitjean, P., Ledoux, C, Vreeswijk, P., Smette, 

A., Meurs, E. J. A. 2009, A&A, 503, 771 
Watson, D., Hjorth, J., Fynbo, J. P. U., Jakobsson, P., Foley, S., 

Sollerman, J., Wijers, R. A. M. J. 2007, ApJ, 660, L101 



a key ingredient for dust production (Draine 2003), our find- 
ings are consistent with a picture in which the darkness of 
a GRB is in most cases due to absorption by circumburst 
material. 
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